Stem cell-derived cardiomyocyte heterogeneity confounds electrophysiological insights.

Human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) offer potential as an in vitro model for studying drug cardiotoxicity and patient-specific cardiovascular disease. The inherent electrophysiological heterogeneity of these cells limits the depth of insights that can be drawn from well-designed experiments. In this review, we provide our perspective on some sources and the consequences of iPSC-CM heterogeneity. We demonstrate the extent of heterogeneity in the literature and explain how such heterogeneity is exacerbated by patch-clamp experimental artifacts in the manual and automated set-up. Finally, we discuss how this heterogeneity, caused by both intrinsic and extrinsic factors, limits our ability to build digital twins of patient-derived cardiomyocytes.

Single-cell ionic current phenotyping explains stem cell-derived cardiomyocyte action potential morphology.

Human induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) are a promising tool to study arrhythmia-related factors, but the variability of action potential (AP) recordings from these cells limits their use as an in vitro model. In this study, we use recently published brief (10 s), dynamic voltage-clamp (VC) data to provide mechanistic insights into the ionic currents contributing to AP heterogeneity; we call this approach rapid ionic current phenotyping (RICP). Features of this VC data were correlated to AP recordings from the same cells, and we used computational models to generate mechanistic insights into cellular heterogeneity. This analysis uncovered several interesting links between AP morphology and ionic current density: both L-type calcium and sodium currents contribute to upstroke velocity, rapid delayed rectifier K+ current is the main determinant of the maximal diastolic potential, and an outward current in the activation range of slow delayed rectifier K+ is the main determinant of AP duration. Our analysis also identified an outward current in several cells at 6 mV that is not reproduced by iPSC-CM mathematical models but contributes to determining AP duration. RICP can be used to explain how cell-to-cell variability in ionic currents gives rise to AP heterogeneity. Because of its brief duration (10 s) and ease of data interpretation, we recommend the use of RICP for single-cell patch-clamp experiments that include the acquisition of APs.NEW & NOTEWORTHY We present rapid ionic current phenotyping (RICP), a current quantification approach based on an optimized voltage-clamp protocol. The method captures a rich snapshot of the ionic current dynamics, providing quantitative information about multiple currents (e.g., ICa,L, IKr) in the same cell. The protocol helped to identify key ionic determinants of cellular action potential heterogeneity in iPSC-CMs. This included unexpected results, such as the critical role of IKr in establishing the maximum diastolic potential.

Optimization of a cardiomyocyte model illuminates role of increased INa,L in repolarization reserve.

Cardiac ion currents may compensate for each other when one is compromised by a congenital or drug-induced defect. Such redundancy contributes to a robust repolarization reserve that can prevent the development of lethal arrhythmias. Most efforts made to describe this phenomenon have quantified contributions by individual ion currents. However, it is important to understand the interplay between all major ion-channel conductances, as repolarization reserve is dependent on the balance between all ion currents in a cardiomyocyte. Here, a genetic algorithm was designed to derive profiles of nine ion-channel conductances that optimize repolarization reserve in a mathematical cardiomyocyte model. Repolarization reserve was quantified using a previously defined metric, repolarization reserve current, i.e., the minimum constant current to prevent normal action potential repolarization in a cell. The optimization improved repolarization reserve current up to 84% compared to baseline in a human adult ventricular myocyte model and increased resistance to arrhythmogenic insult. The optimized conductance profiles were not only characterized by increased repolarizing current conductances but also uncovered a previously unreported behavior by the late sodium current. Simulations demonstrated that upregulated late sodium increased action potential duration, without compromising repolarization reserve current. The finding was generalized to multiple models. Ultimately, this computational approach, in which multiple currents were studied simultaneously, illuminated mechanistic insights into how the metric's magnitude could be increased and allowed for the unexpected role of late sodium to be elucidated.NEW & NOTEWORTHY Genetic algorithms are typically used to fit models or extract desired parameters from data. Here, we use the tool to produce a ventricular cardiomyocyte model with increased repolarization reserve. Since arrhythmia mitigation is dependent on multiple cardiac ion-channel conductances, study using a comprehensive, unbiased, and systems-level approach is important. The use of this optimization strategy allowed us to find robust profiles that illuminated unexpected mechanistic determinants of key ion-channel conductances in repolarization reserve.

Leak current, even with gigaohm seals, can cause misinterpretation of stem cell-derived cardiomyocyte action potential recordings.

AIMS: Human-induced pluripotent stem cell-derived cardiomyocytes (iPSC-CMs) have become an essential tool to study arrhythmia mechanisms. Much of the foundational work on these cells, as well as the computational models built from the resultant data, has overlooked the contribution of seal-leak current on the immature and heterogeneous phenotype that has come to define these cells. The aim of this study is to understand the effect of seal-leak current on recordings of action potential (AP) morphology. METHODS AND RESULTS: Action potentials were recorded in human iPSC-CMs using patch clamp and simulated using previously published mathematical models. Our in silico and in vitro studies demonstrate how seal-leak current depolarizes APs, substantially affecting their morphology, even with seal resistances (Rseal) above 1 GΩ. We show that compensation of this leak current is difficult due to challenges with obtaining accurate measures of Rseal during an experiment. Using simulation, we show that Rseal measures (i) change during an experiment, invalidating the use of pre-rupture values, and (ii) are polluted by the presence of transmembrane currents at every voltage. Finally, we posit that the background sodium current in baseline iPSC-CM models imitates the effects of seal-leak current and is increased to a level that masks the effects of seal-leak current on iPSC-CMs. CONCLUSION: Based on these findings, we make recommendations to improve iPSC-CM AP data acquisition, interpretation, and model-building. Taking these recommendations into account will improve our understanding of iPSC-CM physiology and the descriptive ability of models built from such data.

Rapid ionic current phenotyping (RICP) identifies mechanistic underpinnings of iPSC-CM AP heterogeneity.

As a renewable, easily accessible, human-derived in vitro model, human induced pluripotent stem cell derived cardiomyocytes (iPSC-CMs) are a promising tool for studying arrhythmia-related factors, including cardiotoxicity and congenital proarrhythmia risks. An oft-mentioned limitation of iPSC-CMs is the abundant cell-to-cell variability in recordings of their electrical activity. Here, we develop a new method, rapid ionic current phenotyping (RICP), that utilizes a short (10 s) voltage clamp protocol to quantify cell-to-cell heterogeneity in key ionic currents. We correlate these ionic current dynamics to action potential recordings from the same cells and produce mechanistic insights into cellular heterogeneity. We present evidence that the L-type calcium current is the main determinant of upstroke velocity, rapid delayed rectifier K+ current is the main determinant of the maximal diastolic potential, and an outward current in the excitable range of slow delayed rectifier K+ is the main determinant of action potential duration. We measure an unidentified outward current in several cells at 6 mV that is not recapitulated by iPSC-CM mathematical models but contributes to determining action potential duration. In this way, our study both quantifies cell-to-cell variability in membrane potential and ionic currents, and demonstrates how the ionic current variability gives rise to action potential heterogeneity. Based on these results, we argue that iPSC-CM heterogeneity should not be viewed simply as a problem to be solved but as a model system to understand the mechanistic underpinnings of cellular variability.

An in silico-in vitro pipeline for drug cardiotoxicity screening identifies ionic pro-arrhythmia mechanisms.

BACKGROUND AND PURPOSE: Before advancing to clinical trials, new drugs are screened for their pro-arrhythmic potential using a method that is overly conservative and provides limited mechanistic insight. The shortcomings of this approach can lead to the mis-classification of beneficial drugs as pro-arrhythmic. EXPERIMENTAL APPROACH: An in silico-in vitro pipeline was developed to circumvent these shortcomings. A computational human induced pluripotent stem cell-derived cardiomyocyte (iPSC-CM) model was used as part of a genetic algorithm to design experiments, specifically electrophysiological voltage clamp (VC) protocols, to identify which of several cardiac ion channels were blocked during in vitro drug studies. Such VC data, along with dynamically clamped action potentials (AP), were acquired from iPSC-CMs before and after treatment with a control solution or a low- (verapamil), intermediate- (cisapride or quinine) or high-risk (quinidine) drug. KEY RESULTS: Significant AP prolongation (a pro-arrhythmia marker) was seen in response to quinidine and quinine. The VC protocol identified block of IKr (a source of arrhythmias) by all strong IKr blockers, including cisapride, quinidine and quinine. The protocol also detected block of ICaL by verapamil and Ito by quinidine. Further demonstrating the power of the approach, the VC data uncovered a previously unidentified If block by quinine, which was confirmed with experiments using a HEK-293 expression system and automated patch-clamp. CONCLUSION AND IMPLICATIONS: We developed an in silico-in vitro pipeline that simultaneously identifies pro-arrhythmia risk and mechanism of ion channel-blocking drugs. The approach offers a new tool for evaluating cardiotoxicity during preclinical drug screening.

A dual SHOX2:GFP; MYH6:mCherry knockin hESC reporter line for derivation of human SAN-like cells.

The sinoatrial node (SAN) is the primary pacemaker of the heart. The human SAN is poorly understood due to limited primary tissue access and limitations in robust in vitro derivation methods. We developed a dual SHOX2:GFP; MYH6:mCherry knockin human embryonic stem cell (hESC) reporter line, which allows the identification and purification of SAN-like cells. Using this line, we performed several rounds of chemical screens and developed an efficient strategy to generate and purify hESC-derived SAN-like cells (hESC-SAN). The derived hESC-SAN cells display molecular and electrophysiological characteristics of bona fide nodal cells, which allowed exploration of their transcriptional profile at single-cell level. In sum, our dual reporter system facilitated an effective strategy for deriving human SAN-like cells, which can potentially be used for future disease modeling and drug discovery.

Validation of quantitative measure of repolarization reserve as a novel marker of drug induced proarrhythmia.

Repolarization reserve, the robustness of a cell to repolarize even when one of the repolarization mechanisms is failing, has been described qualitatively in terms of ionic currents, but has not been quantified by a generic metric that is applicable to drug screening. Prolonged repolarization leading to repolarization failure is highly arrhythmogenic. It may lead to ventricular tachycardia caused by triggered activity from early afterdepolarizations (EADs), or it may promote the occurrence of unidirectional conduction block and reentry. Both types of arrhythmia may deteriorate into ventricular fibrillation (VF) and death. We define the Repolarization Reserve Current (RRC) as the minimum constant current necessary to prevent normal repolarization of a cell. After developing and testing RRC for nine computational ionic models of various species, we applied it experimentally to atrial and ventricular human induced pluripotent stem cell-derived cardiomyocyte (hiPSC-CM), and isolated guinea-pig ventricular cardiomyocytes. In simulations, repolarization was all-or-none with a precise, model-dependent critical RRC, resulting in a discrete shift in the Action Potential Duration (APD) - RRC relation, in the occurrence of EADs and repolarization failure. These data were faithfully reproduced in cellular experiments. RRC allows simple, fast, unambiguous quantification of the arrhythmogenic propensity in cardiac cells of various origins and species without the need of prior knowledge of underlying currents and is suitable for high throughput applications, and personalized medicine applications.

Calibration of ionic and cellular cardiac electrophysiology models.

Cardiac electrophysiology models are among the most mature and well-studied mathematical models of biological systems. This maturity is bringing new challenges as models are being used increasingly to make quantitative rather than qualitative predictions. As such, calibrating the parameters within ion current and action potential (AP) models to experimental data sets is a crucial step in constructing a predictive model. This review highlights some of the fundamental concepts in cardiac model calibration and is intended to be readily understood by computational and mathematical modelers working in other fields of biology. We discuss the classic and latest approaches to calibration in the electrophysiology field, at both the ion channel and cellular AP scales. We end with a discussion of the many challenges that work to date has raised and the need for reproducible descriptions of the calibration process to enable models to be recalibrated to new data sets and built upon for new studies. This article is categorized under: Analytical and Computational Methods > Computational Methods Physiology > Mammalian Physiology in Health and Disease Models of Systems Properties and Processes > Cellular Models.

A human embryonic stem cell reporter line for monitoring chemical-induced cardiotoxicity.

AIMS: Human embryonic stem cells (hESCs) can be used to generate scalable numbers of cardiomyocytes (CMs) for studying cardiac biology, disease modelling, drug screens, and potentially for regenerative therapies. A fluorescence-based reporter line will significantly enhance our capacities to visualize the derivation, survival, and function of hESC-derived CMs. Our goal was to develop a reporter cell line for real-time monitoring of live hESC-derived CMs. METHODS AND RESULTS: We used CRISPR/Cas9 to knock a mCherry reporter gene into the MYH6 locus of hESC lines, H1 and H9, enabling real-time monitoring of the generation of CMs. MYH6:mCherry+ cells express atrial or ventricular markers and display a range of cardiomyocyte action potential morphologies. At 20 days of differentiation, MYH6:mCherry+ cells show features characteristic of human CMs and can be used successfully to monitor drug-induced cardiotoxicity and oleic acid-induced cardiac arrhythmia. CONCLUSION: We created two MYH6:mCherry hESC reporter lines and documented the application of these lines for disease modelling relevant to cardiomyocyte biology.

Light-Activated Dynamic Clamp Using iPSC-Derived Cardiomyocytes.

iPSC-derived cardiomyocytes (iPSC-CMs) are a potentially advantageous platform for drug screening because they provide a renewable source of human cardiomyocytes. One obstacle to their implementation is their immature electrophysiology, which reduces relevance to adult arrhythmogenesis. To address this, dynamic clamp is used to inject current representing the insufficient potassium current, IK1, thereby producing more adult-like electrophysiology. However, dynamic clamp requires patch clamp and is therefore low throughput and ill-suited for large-scale drug screening. Here, we use optogenetics to generate such a dynamic-clamp current. The optical dynamic clamp (ODC) uses outward-current-generating opsin, ArchT, to mimic IK1, resulting in more adult-like action potential morphology, similar to IK1 injection via classic dynamic clamp. Furthermore, in the presence of an IKr blocker, ODC revealed expected action potential prolongation and reduced spontaneous excitation. The ODC presented here still requires an electrode to measure Vm but provides a first step toward contactless dynamic clamp, which will not only enable high-throughput screening but may also allow control within multicellular iPSC-CM formats to better recapitulate adult in vivo physiology.

Overexpression of Map3k7 activates sinoatrial node-like differentiation in mouse ES-derived cardiomyocytes.

In vivo, cardiomyocytes comprise a heterogeneous population of contractile cells defined by unique electrophysiologies, molecular markers and morphologies. The mechanisms directing myocardial cells to specific sub-lineages remain poorly understood. Here we report that overexpression of TGFß-Activated Kinase (TAK1/Map3k7) in mouse embryonic stem (ES) cells faithfully directs myocardial differentiation of embryoid body (EB)-derived cardiac cells toward the sinoatrial node (SAN) lineage. Most cardiac cells in Map3k7-overexpressing EBs adopt markers, cellular morphologies, and electrophysiological behaviors characteristic of the SAN. These data, in addition to the fact that Map3k7 is upregulated in the sinus venous-the source of cells for the SAN-suggest that Map3k7 may be an endogenous regulator of the SAN fate.

Slow [Na+]i dynamics impacts arrhythmogenesis and spiral wave reentry in cardiac myocyte ionic model.

Accumulation of intracellular Na+ is gaining recognition as an important regulator of cardiac myocyte electrophysiology. The intracellular Na+ concentration can be an important determinant of the cardiac action potential duration, can modulate the tissue-level conduction of excitation waves, and can alter vulnerability to arrhythmias. Mathematical models of cardiac electrophysiology often incorporate a dynamic intracellular Na+ concentration, which changes much more slowly than the remaining variables. We investigated the dependence of several arrhythmogenesis-related factors on [Na+]i in a mathematical model of the human atrial action potential. In cell simulations, we found that [Na+]i accumulation stabilizes the action potential duration to variations in several conductances and that the slow dynamics of [Na+]i impacts bifurcations to pro-arrhythmic afterdepolarizations, causing intermittency between different rhythms. In long-lasting tissue simulations of spiral wave reentry, [Na+]i becomes spatially heterogeneous with a decreased area around the spiral wave rotation center. This heterogeneous region forms a functional anchor, resulting in diminished meandering of the spiral wave. Our findings suggest that slow, physiological, rate-dependent variations in [Na+]i may play complex roles in cellular and tissue-level cardiac dynamics.

Hard real-time closed-loop electrophysiology with the Real-Time eXperiment Interface (RTXI).

The ability to experimentally perturb biological systems has traditionally been limited to static pre-programmed or operator-controlled protocols. In contrast, real-time control allows dynamic probing of biological systems with perturbations that are computed on-the-fly during experimentation. Real-time control applications for biological research are available; however, these systems are costly and often restrict the flexibility and customization of experimental protocols. The Real-Time eXperiment Interface (RTXI) is an open source software platform for achieving hard real-time data acquisition and closed-loop control in biological experiments while retaining the flexibility needed for experimental settings. RTXI has enabled users to implement complex custom closed-loop protocols in single cell, cell network, animal, and human electrophysiology studies. RTXI is also used as a free and open source, customizable electrophysiology platform in open-loop studies requiring online data acquisition, processing, and visualization. RTXI is easy to install, can be used with an extensive range of external experimentation and data acquisition hardware, and includes standard modules for implementing common electrophysiology protocols.

Global Optimization of Ventricular Myocyte Model to Multi-Variable Objective Improves Predictions of Drug-Induced Torsades de Pointes.

In silico cardiac myocyte models present powerful tools for drug safety testing and for predicting phenotypical consequences of ion channel mutations, but their accuracy is sometimes limited. For example, several models describing human ventricular electrophysiology perform poorly when simulating effects of long QT mutations. Model optimization represents one way of obtaining models with stronger predictive power. Using a recent human ventricular myocyte model, we demonstrate that model optimization to clinical long QT data, in conjunction with physiologically-based bounds on intracellular calcium and sodium concentrations, better constrains model parameters. To determine if the model optimized to congenital long QT data better predicts risk of drug-induced long QT arrhythmogenesis, in particular Torsades de Pointes risk, we tested the optimized model against a database of known arrhythmogenic and non-arrhythmogenic ion channel blockers. When doing so, the optimized model provided an improved risk assessment. In particular, we demonstrate an elimination of false-positive outcomes generated by the baseline model, in which simulations of non-torsadogenic drugs, in particular verapamil, predict action potential prolongation. Our results underscore the importance of currents beyond those directly impacted by a drug block in determining torsadogenic risk. Our study also highlights the need for rich data in cardiac myocyte model optimization and substantiates such optimization as a method to generate models with higher accuracy of predictions of drug-induced cardiotoxicity.

Differential roles of two delayed rectifier potassium currents in regulation of ventricular action potential duration and arrhythmia susceptibility.

KEY POINTS: Arrhythmias result from disruptions to cardiac electrical activity, although the factors that control cellular action potentials are incompletely understood. We combined mathematical modelling with experiments in heart cells from guinea pigs to determine how cellular electrical activity is regulated. A mismatch between modelling predictions and the experimental results allowed us to construct an improved, more predictive mathematical model. The balance between two particular potassium currents dictates how heart cells respond to perturbations and their susceptibility to arrhythmias. ABSTRACT: Imbalances of ionic currents can destabilize the cardiac action potential and potentially trigger lethal cardiac arrhythmias. In the present study, we combined mathematical modelling with information-rich dynamic clamp experiments to determine the regulation of action potential morphology in guinea pig ventricular myocytes. Parameter sensitivity analysis was used to predict how changes in ionic currents alter action potential duration, and these were tested experimentally using dynamic clamp, a technique that allows for multiple perturbations to be tested in each cell. Surprisingly, we found that a leading mathematical model, developed with traditional approaches, systematically underestimated experimental responses to dynamic clamp perturbations. We then re-parameterized the model using a genetic algorithm, which allowed us to estimate ionic current levels in each of the cells studied. This unbiased model adjustment consistently predicted an increase in the rapid delayed rectifier K+ current and a drastic decrease in the slow delayed rectifier K+ current, and this prediction was validated experimentally. Subsequent simulations with the adjusted model generated the clinically relevant prediction that the slow delayed rectifier is better able to stabilize the action potential and suppress pro-arrhythmic events than the rapid delayed rectifier. In summary, iterative coupling of simulations and experiments enabled novel insight into how the balance between cardiac K+ currents influences ventricular arrhythmia susceptibility.

Applications of Dynamic Clamp to Cardiac Arrhythmia Research: Role in Drug Target Discovery and Safety Pharmacology Testing.

Dynamic clamp, a hybrid-computational-experimental technique that has been used to elucidate ionic mechanisms underlying cardiac electrophysiology, is emerging as a promising tool in the discovery of potential anti-arrhythmic targets and in pharmacological safety testing. Through the injection of computationally simulated conductances into isolated cardiomyocytes in a real-time continuous loop, dynamic clamp has greatly expanded the capabilities of patch clamp outside traditional static voltage and current protocols. Recent applications include fine manipulation of injected artificial conductances to identify promising drug targets in the prevention of arrhythmia and the direct testing of model-based hypotheses. Furthermore, dynamic clamp has been used to enhance existing experimental models by addressing their intrinsic limitations, which increased predictive power in identifying pro-arrhythmic pharmacological compounds. Here, we review the recent advances of the dynamic clamp technique in cardiac electrophysiology with a focus on its future role in the development of safety testing and discovery of anti-arrhythmic drugs.

Illuminating Myocyte-Fibroblast Homotypic and Heterotypic Gap Junction Dynamics Using Dynamic Clamp.

Fibroblasts play a significant role in the development of electrical and mechanical dysfunction of the heart; however, the underlying mechanisms are only partially understood. One widely studied mechanism suggests that fibroblasts produce excess extracellular matrix, resulting in collagenous septa that slow propagation, cause zig-zag conduction paths, and decouple cardiomyocytes, resulting in a substrate for cardiac arrhythmia. An emerging mechanism suggests that fibroblasts promote arrhythmogenesis through direct electrical interactions with cardiomyocytes via gap junction (GJ) channels. In the heart, three major connexin (Cx) isoforms, Cx40, Cx43, and Cx45, form GJ channels in cell-type-specific combinations. Because each Cx is characterized by a unique time- and transjunctional voltage-dependent profile, we investigated whether the electrophysiological contributions of fibroblasts would vary with the specific composition of the myocyte-fibroblast (M-F) GJ channel. Due to the challenges of systematically modifying Cxs in vitro, we coupled native cardiomyocytes with in silico fibroblast and GJ channel electrophysiology models using the dynamic-clamp technique. We found that there is a reduction in the early peak of the junctional current during the upstroke of the action potential (AP) due to GJ channel gating. However, effects on the cardiomyocyte AP morphology were similar regardless of the specific type of GJ channel (homotypic Cx43 and Cx45, and heterotypic Cx43/Cx45 and Cx45/Cx43). To illuminate effects at the tissue level, we performed multiscale simulations of M-F coupling. First, we developed a cell-specific model of our dynamic-clamp experiments and investigated changes in the underlying membrane currents during M-F coupling. Second, we performed two-dimensional tissue sheet simulations of cardiac fibrosis and incorporated GJ channels in a cell type-specific manner. We determined that although GJ channel gating reduces junctional current, it does not significantly alter conduction velocity during cardiac fibrosis relative to static GJ coupling. These findings shed more light on the complex electrophysiological interplay between cardiac fibroblasts and myocytes.

Improving cardiomyocyte model fidelity and utility via dynamic electrophysiology protocols and optimization algorithms.

Mathematical models of cardiac electrophysiology are instrumental in determining mechanisms of cardiac arrhythmias. However, the foundation of a realistic multiscale heart model is only as strong as the underlying cell model. While there have been myriad advances in the improvement of cellular-level models, the identification of model parameters, such as ion channel conductances and rate constants, remains a challenging problem. The primary limitations to this process include: (1) such parameters are usually estimated from data recorded using standard electrophysiology voltage-clamp protocols that have not been developed with model building in mind, and (2) model parameters are typically tuned manually to subjectively match a desired output. Over the last decade, methods aimed at overcoming these disadvantages have emerged. These approaches include the use of optimization or fitting tools for parameter estimation and incorporating more extensive data for output matching. Here, we review recent advances in parameter estimation for cardiomyocyte models, focusing on the use of more complex electrophysiology protocols and global search heuristics. We also discuss future applications of such parameter identification, including development of cell-specific and patient-specific mathematical models to investigate arrhythmia mechanisms and predict therapy strategies.

Computational Approaches to Understanding the Role of Fibroblast-Myocyte Interactions in Cardiac Arrhythmogenesis.

The adult heart is composed of a dense network of cardiomyocytes surrounded by nonmyocytes, the most abundant of which are cardiac fibroblasts. Several cardiac diseases, such as myocardial infarction or dilated cardiomyopathy, are associated with an increased density of fibroblasts, that is, fibrosis. Fibroblasts play a significant role in the development of electrical and mechanical dysfunction of the heart; however the underlying mechanisms are only partially understood. One widely studied mechanism suggests that fibroblasts produce excess extracellular matrix, resulting in collagenous septa. These collagenous septa slow propagation, cause zig-zag conduction paths, and decouple cardiomyocytes resulting in a substrate for arrhythmia. Another emerging mechanism suggests that fibroblasts promote arrhythmogenesis through direct electrical interactions with cardiomyocytes via gap junctions. Due to the challenges of investigating fibroblast-myocyte coupling in native cardiac tissue, computational modeling and in vitro experiments have facilitated the investigation into the mechanisms underlying fibroblast-mediated changes in cardiomyocyte action potential morphology, conduction velocity, spontaneous excitability, and vulnerability to reentry. In this paper, we summarize the major findings of the existing computational studies investigating the implications of fibroblast-myocyte interactions in the normal and diseased heart. We then present investigations from our group into the potential role of voltage-dependent gap junctions in fibroblast-myocyte interactions.